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Abstract: Distinguishing 7 -ray and hadron initiated extensive air showers (EAS) based on lateral distribution 
of electrons has been studied by detailed Monte Carlo (MC) simulations. The possibility of using the local age 
parameter (LAP) of EAS for the gamma-hadron separation has been explored. It is found that separating 7 -ray 
and hadron induced EAS on the basis of LAP can be useful for surface detector experiments those have no 
reliable muon measurement facilities. 
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1 Introduction 

The important features of the primary energy spectrum vis- 
a-vis the chemical composition of the cosmic rays (CRs) at 
ultra high energy (UHE) and extremely high energy can be 
understood from a complete study of CR air-shower. The 
direct measurements of primary CR flux nearly 100 TeV 
and above are impractical because of very low and sharply 
falling flux, but has to be inferred from observations of 
extensive air showers. 

To extract information about primary CRs from the mea¬ 
surements of various experiments still require detailed M- 
C simulations of the shower development as a basis of the 
data analysis and interpretation. The MC simulations con¬ 
sider the evolution of EAS in the atmosphere initiated by 
different energetic particles. 

In the field of CR air shower physics, the discrimina¬ 
tion of 7 -ray induced air showers from hadron-induced air 
showers is a challenging problem that still needs more at¬ 
tention m. The high energy 7 -rays are considered to be 
the by-product of hadronic components of primary CRs: a 
fraction of accelerated hadronic cosmic rays are likely to 
produce 7 -rays through A-resonances, bremsstrahlung, in¬ 
verse Compton etc. processes at or very close to the source 
site by interacting with the ambient matter El- The 7 -ray 
detection suffers from the huge background constituted 
by ordinary CRs (hadrons) in the GeV-TeV energy region. 
To observe astrophysical sources (emitting undeflected 7 - 
rays and assumed as point-like objects) and to study the 
anisotropy properties of primary CRs, one has to eliminate 
isotropically distributed CRs. 

Usually poor muon content is considered as the signa¬ 
ture of 7 -ray initiated EAS. In order to select a 7 -ray show¬ 
er based on such a criterion, an EAS array needs to be 
equipped with muon detectors covering very large area 
which is economically very challenging and such a facility 
is rarely available. So one has to look for some other pri¬ 
mary mass sensitive observables based on which 7 -ray ini¬ 
tiated EASs can be separated out without the need of large 
area muon detector. The lateral shower age, which is essen¬ 
tially the slope of the lateral density distribution of elec¬ 
trons in EAS reflects the developmental stage of EAS and 
hence it can be used as distinguishing parameter. 


The selection of 7 -ray initiated showers is attempted 
here by employing two different approaches - the Method 

1 and the Method II. In the former approach, we have taken 
single (r-independent) lateral age parameter i.e. ,v as the 
7 -ray separation parameter. Experimentally it is observed 
that the NKG function with a single lateral age is insuf¬ 
ficient to describe the lateral density distribution of EAS 
electrons properly at all distances, which implies that the 
lateral age changes with the radial distance. Subsequently, 
the notion of local lateral shower age parameter (LAP) was 
introduced (3[ which is in essence the lateral age at a point. 
Since experimental electron density data in EAS may fluc¬ 
tuate considerably at a particular radial distance, instead of 
taking LAP at any particular point we take an average LAP 
between 50 m (minimum) and 300 m (maximum) and in 
Method II, we employed this average LAP (s; oca /) to select 
7 -ray showers. 

The present study has been performed mainly at the 
geographical location of ARGO-YBJ (4) (latitude 30.11° 
N, longitude 90.53° E, 4300 m a.s.l.). This is because 
the experiment offers a full coverage array and hence can 
measure radial density distribution of electrons in EAS 
with great accuracy, which in turn provides an opportunity 
to estimate LAP accurately. However, ARGO-YBJ has not 
yet studied the radial variation of the lateral shower age 
while a few other experiments, such as Akeno 0 and 
NBU 0, successfully tested the predicted radial variation 
of LAP. We simulated a few events at the geographical 
location of Akeno and compared with the observations 
to demonstrate again the importance of considering LAP 
instead of single shower age. To ensure that the conclusion 
of the present work is robust, we consider two high energy 
interaction models, QGSJet 01 v.lc Q and EPOS - 1.99 
QD and found that the present findings do not have any 
strong dependence on the choice of interaction model. 

2 Shower age parameters from the cascade 
theory 

Nishimura and Kamata J9l solved numerically the 3- 
dimensional shower equations in Approximation B (elec¬ 
trons suffer constant amount of collision loss in a radiation 
length which is equal to the critical energy) to obtain the 
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lateral distribution of electrons propagating in a medium 
of constant density. The results obtained on lateral densi¬ 
ty distribution of cascade particles by Nishimura and Ka- 
mata can be approximated by the well known Nishimura- 
Kamata-Greisen (NKG) structure function proposed by 
Greisen flOl . given by 


f(r) = C{s 1 _){r/r m ) s ± 2 (1 + r/r m ) s ± 45 , (1) 

where the normalization factor C(.v |). 

The normalization of f[r) implies for the electron densi¬ 
ty p ( r ) = N e f(r). But such a widely used relation does not 
hold if .v varies with r, as noted in Akeno (5), NBU |6| 
and some other observations 0. Consequent upon, the lat¬ 
eral shower age is found to vary with radial distance exper¬ 
imentally. To handle the situation a method was developed 
by Capdevielle et al introducing the notion of local age pa¬ 
rameter (LAP) 0- From two neighboring points, i and j, 
we can give a LAP for any distribution f(x) (where x = j-) 
which characterizes the best fit by a NKG-type function in 
[xuxj] : 


Slocal (}■>./) 
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Here, F tj = /(r,-)//(r ; -), X^r^rj, and Y l j=(x i +])/(xj+]). 
More generally, if r,- -a rj, this suggests the definition of 
the LAP si oca i (x) (or s /oca /(r)) at each point: 
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Function fNKc( r ) with s=s/ oca i(r) can be used to fit f{r) 
in the neighborhood of r. 

The identification si oca ,{r) = si oca i(i.j) for r = 
remains valid for the experimental distributions (taking 
Fjj = p(r,)/p(rj)) as far as they are approximated by 
monotonic decreasing functions versus distance. 

The behavior of LAP on experimental lateral distribu¬ 
tions was found in accordance with the prediction lUTl 
which was reaffirmed by the Akeno observations 0. The 
stated method was validated by the rapporteurs of the ICR- 
C from 1981 to 1985 lfl2l. 

3 Simulation of EAS 

The EAS events are simulated by coupling the high energy 
(above 80 GeV/n) hadronic interaction models QGSJet 01 
version lc 0 and EPOS 1.99 0, and the low energy (be¬ 
low 80 GeV/n) hadronic interaction model GHEISHA (ver¬ 
sion 2002d) ifBl in the framework of the CORSIKA Monte 
Carlo program version 6.970 M . For the electromagnetic 
part the EGS4 lfl5ll program library has been used. 

The MC simulation data library consists of considerable 
amounts of EAS events each for p, Fe and y-ray at the 
ARGO-YBJ and Akeno levels. A mixed sample has been 
prepared from the generated showers taking 37% p, 37% 
Fe, 26% y-ray events for better understanding of EAS ob¬ 
servational results. 

4 Estimation of lateral shower age 

The simulated electron density data have been analyzed 
in two different methods to obtain shower age parameters 


Fig. 1: [a] The radial variation of LAP obtained at the 

ARGO-YBJ altitude for p, Fe and y-ray with hadronic 
interaction models QGSJet and EPOS, [b] Same as figure 
la but at the Akeno level with QGSJet model, compared 
with experimental data. The solid and dashed lines parallel 
to the x-axis indicate ,v_ L for Fe and p. 

s_l and si oca i {r). First, following the traditional approach 
we estimate .y by fitting the density data with the NKG 
structure function. The error in estimating ,v | in the N e 
range 1.5 x 10 4 — 4.25 x 10 5 has been found as ±0.03. 

Secondly, exploiting Eq.(2) we directly estimated LAP 
for each individual event. In this analysis, the error of the 
LAP for EAS with the primary energy in the PeV range 
remains within 0.05 for 10 m < r < 250 m, whereas for 
r < 10 m or when r > 300 m the error of the LAP is found 
to be higher, about 0.1. 

It has been shown recently that LAP initially decreases 
with radial distance, reaches a minimum around 50 m, 
then starts increasing with radial distance, attains a local 
maximum around 300 m and decreases again thereafter 
mu. We observed such a characteristic variation of LAP 
with radial distance at ARGO-YBJ location as shown in 
the top figure Q] for proton, iron and y-ray primaries for 
both QGSJet and EPOS models. 

As already mentioned, Akeno group studied the radial 
variation of LAP experimentally 0- We compared our 
simulation results with Akeno observations in the bottom 
figure [T| for proton and Fe primaries. The single lateral 
shower age for proton and iron primaries are also given 
(solid and dashed lines parallel to the x-axis) in the figure 
for comparison. 

Since air shower measurements are subjected to large 
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Fig. 2: [a] Variation of mean lateral shower age with show¬ 
er size at ARGO-YBJ level; [b] Variation of mean local 
age parameter with shower size at ARGO-YBJ level. 

fluctuations, instead of LAP at a particular radial distance 
we consider for each event a mean LAP (< si oca i >), which 
is the average of LAPs for several small distance bands 
(n,fj) over the radial distance between 50 m to 300 m. 
For the purpose of averaging, distance bands are taken in 
constant steps on the logarithmic scale. The radial distance 
band from 50 m to 300 m is chosen because the positions 
of local minimum and maximum at 50 m and 300 m are 
nearly universal, independently of primary energy lfl6l . 

5 Potentiality of Method II over Method I 

The inadequacy of a single (constant) lateral age parameter 
to describe the experimental lateral distribution of EAS 
electrons properly at all distances has been noted in several 
experimental observations. So when a single constant age 
is assigned to an EAS event, the discriminating power of 
that parameter on primary masses somewhat becomes dull; 
it still can distinguish iron initiated showers from proton 
induced showers but can’t effectively separate out y-ray 
showers from EASs generated by primary protons and this 
is clear from the comparison between the top and bottom 
figures in [2] 

6 Selection of y-ray component using 
LAP:Method II 

In the primary CR flux, the percentage of y-ray flux is very 
very small, of the order of 0.001%. To separate such a s¬ 


Fig. 3: Distribution of simulated events in mixture-II in a 
primary energy range based on s/ oea / and N e at ARGO- 
YBJ.Dependence on interaction models is also shown. 

mall fraction of y-ray component from primary CRs in a 
real experiment employing method II, it would be nice if 
we could prepare a mixture with such a small percent of 
y-rays and see whether the y-ray events can be extracted 
out or not. But due to limited statistics we could not do 
that. Instead, we made three sub-mixtures of type II main¬ 
taining the ratio of primaries as 37% p, 37% Fe and 26% 
y-ray for three different combinations of zenith angle, pri¬ 
mary energy and shower size range (scatter plots in figures 
[T|i. These exercises are done for both the high energy inter¬ 
action models QGSJet and EPOS. From each sub-mixture, 
we tried to separate y-ray showers out exploiting si oca i cor¬ 
responding to different /V,,. 

For a small enough cut value of si oca i we have found 
poor acceptance for y-ray induced EAS and very good re¬ 
jection of background, whereas for a large enough value of 
siocai the situation is found to be completely reversed. Cut 
values of s/ ora / lie between the two extremes offered differ¬ 
ent acceptances and rejections of y-ray and background re¬ 
spectively. Similarly a proper cut on N e is also needed for 
selection of y-ray induced EAS and the rejection of back¬ 
ground. In astronomical signal selection, optimal cut to s- 
election parameters (maximizing the y-ray efficiency and 
minimizing the background contamination) is usually set 
by numerical maximization of the quality factor Q defined 
by 


V £bk g V~ £ bkg ) 

where Ey and respectively denote the acceptances 
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Siocai 

1.05 

1.08 

1.10 

1.12 

1.14 

1.15 

1.16 

1.18 

1.20 

e 7 

0.146 

0.237 

0.348 

0.510 

0.657 

0.707 

0.778 

0.864 

0.934 

^bkg 

0.007 

0.022 

0.029 

0.036 

0.058 

0.065 

0.109 

0.145 

0.254 

Q 

1.74 

1.61 

2.05 

2.69 

2.73 

2.77 

2.36 

2.27 

1.86 


Table 1: The quality factor at various si oca i cuts using QGSJet model. The primary energy, zenith angle and shower size 
intervals are 200 — 500 TeV, 20° — 30° and (1.5 — 4.0) x 10 5 respectively. 


Model 

E (TeV) 

e (deg.) 

N e X 10 S 

Siocai 

e 7 

^bkg 

Q 

EPOS 

100-200 

5-15 

0 . 6 - 2.0 

1.14 

0.691 

0.072 

2.57 

EPOS 

200-500 

20-30 

O 

1 

i n 

1.15 

0.734 

0.048 

3.35 

QGSJet 

200-500 

20-30 

p 

'xf- 

1 

i n 

1.15 

0.707 

0.065 

2.77 


Table 2: The signal selection parameters at optimal conditions. 


of y-ray and background from the sample using a cut val¬ 
ue of siocai but stands for the rejection of background. 
Q essentially quantifies the gain of significance achieved 
by the separation algorithm. As for an example the quality 
factor is evaluated over a shower size and zenith angle bins 
by varying si oca i which is shown in table (Table Q}. In table 
(Table 0, we have given a chart for quality factors estimat¬ 
ed at optimal conditions in three different situations. 

7 Summary & Conclusions 

In this work, we first explored the r-independent shower 
age as the distinguishing parameter. Later, we attempted 
to separate out y-ray initiated EAS on the basis of mean 
LAP. It is found from the simulation results that y-ray 
induced EASs are younger in terms of the mean LAP and 
hence this parameter (mean LAP) can effectively separate 
out y-ray induced showers from hadronic EAS unlike the 
case of single lateral shower age. Lateral distribution of 
electrons in EAS exhibits universal (primary energy and 
mass independent) behavior in terms of LAP ED- 

An important question is the experimental realization of 
the adopted technique involving LAP. The uncertainty in 
estimating LAP is usually large in normal circumstances 
in comparison to that in lateral shower age as the LAP 
depends on the logarithmic derivative of the density versus 
radial distance. These uncertainties should be small for a 
closely packed air shower array like GRAPES-III at Ooty 
lfl7l or for a full coverage EAS array like ARGO-YBJ |4| . 
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